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Abstract 

We have investigated the behavior of bistable cells made up of four quantum 
dots and occupied by two electrons, in the presence of realistic confinement 
potentials produced by depletion gates on top of a GaAs/AlGaAs heterostruc- 
ture. Such a cell represents the basic building block for logic architectures 
based on the concept of Quantum Cellular Automata (QCA) and of ground 
state computation, which have been proposed as an alternative to traditional 
transistor-based logic circuits. We have focused on the robustness of the op- 
eration of such cells with respect to asymmetries deriving from fabrication 
tolerances. We have developed a 2-D model for the calculation of the electron 
density in a driven cell in response to the polarization state of a driver cell. 
Our method is based on the one-shot Configuration-Interaction technique, 
adapted from molecular chemistry. From the results of our simulations, we 
conclude that an implementation of QCA logic based on simple "hole-arrays" 
is not feasible, because of the extreme sensitivity to fabrication tolerances. As 

1 



an alternative, we propose cells defined by multiple gates, where geometrical 
asymmetries can be compensated for by adjusting the bias voltages. Even 
though not immediately applicable to the implementation of logic gates and 
not suitable for large scale integration, the proposed cell layout should allow 
an experimental demonstration of a chain of QCA cells. 
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I. INTRODUCTION 



Several proposals for the implementation of logic functions and data processing based 
on the concept of Quantum Cellular Automata (QCA) and ground state computation have 
appeared in the literature recently. Tougaw et al. [0] devised an architecture (commonly 
known as "Notre Dame architecture") based on bistable cells which couple electrostatically 
to their nearest neighbors. Each cell consists of four (or five) quantum dots and contains a 
total of two electrons. In the absence of external electric fields, the electrons occupy each 
dot with equal probability. In the presence of a nearby (driver) cell with the two electrons 
forced in the dots along one of the two diagonals, alignment along the same diagonal will 
occur in the driver cell, in the hypothesis of potential barriers large enough as to localize 
the electrons. Based on this principle, it is possible to conceive of bistable cell arrays, in 
which the polarization state enforced at the inputs, at the edges of the arrays, propagates 
in a "domino" fashion || , until the ground state is reached throughout the system, and the 
results of the computation are available in the form of the polarization state of the output 
cells, also located at the edges of the arrays. 

Many issues must be resolved before this computational paradigm can be implemented 
in practice: non-invasive detectors are needed to probe the polarization state of the outputs 
without perturbing the ground state of the system; provision must be made for a time 
evolution of the system that is both fast and reliable; design solutions for the basic cell must 
be developed, yielding a reasonable robustness to fabrication tolerances and compatible with 
large-scale integration on a single chip. 

The focus of the work we are presenting is specifically on the robustness of a single cell, 
coupled to a driver cell, to fabrication tolerances and to asymmetries caused by fluctuations 
in the bias voltages applied to the electrodes defining the cell. In particular, we have studied 
the effect of geometrical and electrical asymmetries on the behavior of QCA cells defined by 
means of lateral electrostatic confinement in a GaAs/AlGaAs heterostructure. 

Although often neglected, robustness to fabrication tolerances and manufacturability 
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are central problems affecting all proposed nanoelectronic devices ||, and their solution is 
a prerequisite for any successful new technology. 

We have considered a basic QCA cell with four quantum dots defined by realistic two- 
dimensional confinement potentials, which are computed from the shape of the metal gates 
at the surface of the heterostructure and the voltages applied to them. Calculation of 
the electron density in such a 2-D artificial molecule is a challenging task. Iterative self- 
consistent methods fail to converge, due to the relatively large electrostatic interaction and 
to the particular symmetries associated with the problem. For this reason, we have de- 
veloped a non-iterative technique based on the Configuration-Interaction method used in 
molecular chemistry JFj . The main drawback of this technique is that it requires rather large 
computational resources, and, thus, we are presenting numerical results only for the case of 
occupancy of two electrons per cell. (Work is currently in progress for the inclusion of up 
to six electrons per cell.) 

With respect to the approaches existing in the literature ||, our method allows a direct 
quantitative estimate of the effects of fabrication and bias tolerances on cell operation and 
does not require the introduction of phenomenological parameters such as the tunneling 
energy, which may be hard to evaluate with a realistic potential. 

In the next section, we provide a detailed statement of the problem we intend to solve 
and describe the cell model together with the technique we have used for the computation of 
the 2D confinement potential. In the third section, the solution of the many-body problem 
is discussed and the one-shot Configuration-Interaction method is introduced. Numerical 
results for a cell occupancy of two electrons and various types of asymmetries are presented 
in Sec. IV, where cell design criteria are also established. 

II. STATEMENT OF THE PROBLEM AND CELL MODEL 

Our aim is to investigate the behavior of two coupled QCA cells, each of which is formed 
by four quantum dots and contains two electrons, as illustrated in Fig. |I| Tunneling between 
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the dots of the same cell is allowed, but not between dots belonging to different cells. If the 
barriers separating the dots within a cell are opaque enough, the electron wave functions will 
be localized, and we shall observe a quasi-classical behavior: the two electrons will repel each 
other and localize in two dots along a diagonal, so as to minimize the electrostatic energy. In 
the case of an isolated, symmetric cell, alignment along either diagonal will occur with equal 
probability. If another (driver) cell is placed in proximity to the cell we are investigating 
(driven cell), as in the case represented in Fig. [I], and the electrons in the driver cell are 
taken to be aligned along a given diagonal, their electric field will destroy the symmetry of 
the driven cell, lifting the degeneracy between the two configurations. This will result in the 
electrons in the driven cell lining up parallel to the electrons in the driver cell. 
Following Lent et al. |2| we define a cell polarization P as 

(gi + Qz) - (ga + Q 4 ) m 
Q1 + Q2 + Q3 + Q4 ' U 

where Qi is the integral of the electron density p(f) over the i— th quadrant of a cell. We 
divide each cell into four quadrants (see Fig. ^) and number them counterclockwise. The 
denominator of Eq. (|Ij) is the total number of electrons in the cell and therefore a constant. 
The procedure for calculating the electron density from the many-electron wave functions 
will be discussed in the following section. 

If the two electrons are aligned along the diagonal corresponding to the first and third 
quadrant, P = 1, while if they are aligned along the other diagonal, P = — 1. For finite 
height barriers, values of P intermediate between -1 and 1 are also possible. We define the 
cell-to-cell response function as the function relating the polarization of the driven cell to 
that of the driver cell. In order to compute the polarization of the driven cell in response to 
each value of the polarization of the driver cell, we need to solve for the electronic structure 
in the driven cell in the presence of the electrostatic potential due to the driver cell. 

For the driven cell we consider a two-dimensional model in the effective mass approxi- 
mation. Such a model is valid as long as the thickness of the dots in the vertical direction, 
corresponding to the thickness of the 2-dimensional electron gas (2DEG) from which they 
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are obtained by lateral confinement, is small compared to their other dimensions. The 2DEG 
is obtained by modulation doping next to a GaAs/AlGaAs heterointerface. 

The two-dimensional confinement potential in the plane of the 2DEG is obtained as a 
result of the action of the metal gates, following the method proposed by Davies et al. 



PJTD], without including, to keep the problem manageable from a computational point of 
view, the self-consistent rearrangement of mobile charge within the heterostructure, except 
for that of the two electrons confined in the cell. In other words, the potential due to the 
gates is evaluated with the analytical expressions described below and is used as the bare 
confinement potential for the definition of a many-body Hamiltonian, whose ground state is 
then evaluated with the Configuration- Interaction method. The occupancy of the cell is fixed 
and corresponds to two electrons for all the numerical results we shall be presenting. We 
assume Fermi level pinning at the surface of the semiconductor, so that the electron-electron 



interaction can be treated by the method of images IfLl j, without requiring the solution of 
the Poisson equation. 

We have considered two basic gate configurations for the definition of the four quantum 
dots that make up a cell. The first configuration, represented in Fig. |3], is rather simple and 
straightforward: the four dots are a consequence of four circular holes in a depleting gate 



covering the surface of the heterostructure ||12|| . The second configuration we have studied 



is described in Fig. [| and is more complex: the four dots are defined by a set of seven metal 
gates which create four minima in the 2-D potential at the 2DEG level. 

The contribution to the bare confinement potential from each gate is computed following 
the method developed in Refs. and |10|. The starting point is the well known result that 



gives the potential inside the semiconductor in terms of its boundary values at the plane of 
the surface: 

xr / v 1 f \z\V g (x',y',0) , , 

V g (x,y,z) = — / dx dy , 2 

2% Js (x — x'Y + (y — y 1 ) 1 + z l 

where z is the vertical coordinate, orthogonal to the heterostructure layers, and the integra- 
tion is performed over the gated surface 5'. Given the applied voltages and the shapes of 



the gates, the confining potential can be easily computed. For simple shapes one can derive 
more compact expressions, by performing some of the integrals in Eq. (0) analytically. For 
the cases considered here, we have used Eq. (3.17) of U for gates with circular holes, and 
the equations in Sect. Ill and Sect. IV of [K| for polygonal gates. 

An example of the results obtained by this procedure is reported in Fig. ^ we show the 
confinement potential at a depth of 50 nm, produced by four holes with a diameter of 90 nm, 
with a distance between the centers of 110 nm. The gate voltage has been set at —0.5 V, in 
order to obtain interdot barriers of reasonable height. 

Following Tougaw et al. || , a uniformly distributed positive background charge has been 
added to each cell. Such a charge does not alter the electrostatic energy splitting between 
the two cell configurations and plays a role, from the point of view of the cell-to-cell response 
function, only if the two cells are very close to each other. In this case, it helps prevent 
alterations of the ground state of the driven cell due to the combined electrostatic repulsion 
of the two electrons in the driver cell, which would tend to push both electrons into the two 
dots on the side further from the driver cell itself. In our model, this positive background 
would not actually be needed to achieve charge neutrality ||, which is already ensured by 
the presence of the gates and of positive charges in the donor layer. 



III. CONFIGURATION-INTERACTION METHOD 

As already mentioned, the solution of a many-electron problem in a potential such as that 
present in a QCA cell is rather challenging. Approaches that are typically used for the sim- 
ulation of quantum dots, based on mean-field approximations of the potential seen by each 



electron and on iterative procedures [|13|-|T7|1, fail to converge when applied to the four-dot 
cell. We can understand the reason for this failure considering that convergence of the self- 
consistent iterative procedures in this class of problems is more and more difficult to achieve 
as the electrostatic interaction increases |18| , and in the presence of quasi-degenerate states. 
As long as the electrostatic interaction is small compared to the confinement energy, it is just 
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a perturbation of the latter, and iterative self-consistent procedures converge monotonically 
to the solution. Otherwise, underrelaxation techniques need to be used, but they also often 
fail when closely spaced states are present, due to symmetries or quasi-symmetries in the 
potential landscape: in this situation the charge will bounce back and forth between quasi- 
degenerate states in consecutive iterations, and convergence is never achieved. Since our 
calculations are currently performed at zero temperature, in order to compare with experi- 
mental results that are typically obtained in the tens of millikelvins range, techniques such 
as the Newton method, which are rather successful in the solution of coupled Schrodinger- 
Poisson problems at finite temperature, cannot be successfully used, due to the sharpness 
of the Fermi function at low temperatures. 

The technique we have implemented is based upon an approach often used in molecu- 
lar chemistry 0, the Configuration-Interaction (CI) method. It consists in approximating 
the iV-electron wave function by a finite linear combination of Slater determinants. Most 
importantly, the CI method is a one-shot method, i.e. it does not involve an iterative cal- 
culation of the wave functions, and, hence, does not suffer from the previously described 
convergence / oscillation problems. 

In order to illustrate the CI method, let us consider an N-electron non-relativistic Hamil- 
tonian with a generic two-body interaction g(fi,fj): 

H = H x + H 2 , (3) 

N *2 N 

Hi = £(-^- V? + V(fi)) = £ h(f t ) , (4) 

2 = 1 Z 171 1=1 

H 2 = £ gift, fj) , with g(fi, fj) = g(fj, n) , (5) 

i < j 

where h is the reduced Planck constant and m is the effective mass of the electron (we 
consider the case of Gallium Arsenide, with m = 0.067mo, mo being the free electron mass). 

As we are concerned with confined systems, we can consider a numerable complete basis 
{<£>i(q)}, where q = (r, s) includes both spatial and spin coordinates, over which the single- 
particle wave function can be expanded . In the following, we shall refer to the ^'s as 
spin-orbitals. Using this basis, we build all the possible independent Slater determinants: 
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"JVfc iqij 

n Nk 



(6) 



W fc (q7v) Vn 2t (qw) ••• <Pn Nk (<lN) 
where the index k labels the Slater determinants and the integer rij k specifies which spin- 
orbital appears in the j— th column of the k— th Slater determinant. 

The infinite set is a complete ortho normal basis for the iV— electron eigenfunctions 
of the Hamiltonian in Eq. (|3|) fll9| ; the i-th eigenfunctions can therefore be written as 



fc=l 

To find the eigenfunctions of H, we must solve the secular equation 

where the infinite-dimensional "Hamiltonian matrix" is 

r N 
H kk , = <$ fc |if |$ fc >> = J n(Hi + H 2 )$ k , J] , 



(7) 



(8) 



(9) 



i=l 



Ei is the i-th eigenvalue of H, the vector c« is made up of the coefficients cn~ and / dqi 
stands for integration over the i-th spatial coordinate and summation over spin orientations. 

In practice, this approach cannot be implemented "exactly" (i.e. choosing a complete, 
infinite set of orthonormal spin-orbitals). In our numerical work, we consider a finite set of 
M spin-orbitals {^(q)}, with j — 1 . . . M. 

Given iV-electrons and M spin-orbitals (with M > N), it is possible to build N"sd 
different Slater determinants, where 



SD 



(10) 



With this choice, the secular equation (|j) becomes an Nsd x Nsd Hermitian eigenvalue 
problem. 



The number of nonzero matrix elements is less than (N"sd) 2 , because all the matrix 
elements between determinants differing by more than two spin-orbitals do vanish, as a 
consequence of the selection rules (Slater's rules 0) presented in Appendix A, where the 
general expressions for the elements of TC are also provided. 

The total number M. of nonzero matrix elements is given by the following expression: 

I A/f — at\ I at\ ( m _ jy \ ( 



M 



NZ 



M-N 
1 



v 1 / 



V 



/ 



v 2 / 



x 



1 + (M — N)N + 



(M - N)(M - N - 1)N(N - 1) 



(11) 



Once the eigenvectors of Eq. @ have been obtained, the A^-electron wave function 
can be computed from Eq. (|?p, and the corresponding electron density is simply given by: 



N J2 |^i(fi,si,q 2 , . . . ,qAr)| 2 rfq 2 • • .rfqAr , 

si J 



(12) 



where si represents the spin orientation coordinate. 

Since we are going to present numerical results for a cell occupancy of two electrons, we 
now focus our attention on the two-electron case. The Hamiltonian for the structure under 
study can be written: 

^ = -^ V ?-^ V 2 + Wn) + Wr 2 ) + V driy (ri) + V dlW (r 2 ) + g(n, f 2 ) (13) 

—,'11 —illi 

where Vcon is the confinement potential computed as described in the previous section, 
V^ r j v is the Coulomb potential due to the charge distribution in the neighboring driver cell, 
and g (7*1,7*2) is the two-body interaction. The two-body interaction includes the effects of 
image charges and is given by: 



1 



1 e 2 



9{ri,r 2 ) 



(14) 



4vre |fi - f 2 \ Ane _ f 2 |2 + ^ z ) 2 2z 

where e = e r eo (e r and eo being the relative permittivity of Gallium Arsenide and the vacuum 
permittivity, respectively), and e is the electron charge. This expression has been obtained 
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by taking one half of the electrostatic energy of the system made up of the electrons and 
their image counterparts, since the energy stored in the image space is purely fictitious. The 



last term of Eq. (HP, due to the interaction of each electron with its own image, yields a 
constant shift of the energy spectrum. 

To apply the CI method to the system described by the Hamiltonian in Eq. flT3|), we start 
by choosing a set of n wave functions {^(r*)}. We shall refer to the {^j(r)}'s as orbitals, 
to distinguish them from the spin-orbitals {(fi(q)}. By combining each of the ipiS with one 
of the two possible spin eigenstates corresponding to the two spin orientations along the 
z— axis, we obtain the set {^(q)} of M = In spin-orbitals. 

With M = 2n spin-orbitals and two electrons, we can construct n(2n — 1) independent 
Slater determinants (see Eq. (|T0|)). In the expansion of Eq. (0), instead, we take into account 
only n 2 Slater determinants, i.e., those composed of spin-orbitals with opposite spins, that 
correspond to states with zero total spin component S z along the z— axis. These states are, 
in general, linear combinations of the singlet state IS" = 0, S z = > and of the triplet state 
\S = l,S z = >. Since we are dealing with a spin- independent Hamiltonian, the triplet 
states IS" = 1, S z = > , \S = 1, S z = ±1 > are degenerate in energy; in addition, the 
corresponding wave functions have the same spatial part. Therefore, no information about 
the energy eigenstates of the system is lost, if only one of the triplet states is used to expand 
the wave function of Eq. (|7|) . 

With the above mentioned restriction on the number of Slater determinants, the matrix 
TC is an n 2 x n 2 Hermitian matrix. We can choose the orbitals ^ to be real, thus making 
the matrix 7i real and symmetric. Finally we notice that 7i is a full matrix, since the $fc's, 
being 2x2 determinants, cannot differ by more than two spin-orbitals. 

The choice of the set of orbitals {ipi} is of crucial importance, since the number n of 
orbitals required to get a satisfactory approximation of the ground state energy and the 
corresponding wave function depends on it. We have used as orbitals the single-electron 
eigenstates for an isolated cell and we have found that with this basis 12 orbitals (i.e. 24 
spin orbitals) are sufficient to get a good accuracy in the results for cell sizes around 200 nm. 
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In order to check the validity of the approximation, we have also performed calculations 
using 24 orbitals, finding that the results are practically identical to those obtained with 
the smaller basis. For larger cell sizes, a larger number of orbitals would be necessary, 
because the electrostatic interaction would grow in importance compared to the confinement 
energy. Therefore, the ground state wave function would deviate further from a single Slater 
determinant built with the single-electron orbitals, and would thus need to be expanded on 
a basis of Slater determinants built from a larger set of single-electron wave functions. 

Let us briefly comment on the relationship of the solution obtained with this method to 
those obtained with other commonly used approaches for the self-consistent solution of the 



electronic structure of quantum dots. In the Hartree method and in the Local Density 
Functional Approximation (LDA) [jl~6|l , the exchange term is ignored or treated in an ap- 
proximate fashion, while the wave functions are obtained at each iteration as eigenfunctions 
of a modified Hamiltonian. When convergence is reached, the obtained wave functions are 
the set of one-electron wave functions minimizing the expectation value of the approximate 
Hamiltonian. 



In the case of the Hartree-Fock method ptl the many-electron wave function is repre 



sented by a single Slater determinant, and when self-consistency is attained, the resulting 
Slater determinant minimizes the expectation value of the Hamiltonian, as in the case of 
the Hartree and LDA methods, but with a properly anti-symmetrized wave function. 

With the Configuration-Interaction method, the wave function is expanded over a basis 
of Slater determinants, which guarantees proper antisymmetrization. The basis functions 
are fixed and the unknowns are the coefficients of the expansion. If the basis were complete 
(comprising an infinite number of spin-orbitals), the solution would be exact and correspond- 
ing to that obtainable by diagonalizing the many-body Hamiltonian. In order to make the 
problem computationally feasible, we must limit the number of basis functions, introducing, 
as a result, some approximation. The difference between the solutions obtained with the 
application of the Hartree-Fock method and of the CI method can thus be summerized as 
follows: with the former we get a single, optimized Slater determinant, while with the latter 
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we obtain an expansion of the solution over a finite basis of Slater determinants which have 
been chosen "a priori". In the presence of strong electron-electron correlation, the wave 
function obtained with the CI method with the inclusion of a reasonable number of basis 
functions is expected to be much closer to the exact solution than the optimized Slater 
determinant resulting from the Hartree-Fock method. 

IV. RESULTS AND DISCUSSION 
A. Single gate configuration with four circular holes 

We start the presentation of the numerical results with the cell-to-cell response function 
obtained for a cell defined by a gate with four circular holes, for a choice of 4, 12, and 24 
basis orbitals (Fig. [|). Each cell is defined by holes with a diameter of 90 nm, placed at the 
corners of a square, with a distance of 110 nm between the hole centers. The applied gate 
bias is —0.5 V and the separation between the centers of the two cells is 400 nm. We found 
that there is no significant difference between the results obtained with the three choices of 
basis elements. In some other cases we have noticed a small difference between the response 
function calculated with 4 basis orbitals and those for 12 and 24 orbitals, which were instead 
found to be practically identical. All the results presented in the following have therefore 
been obtained with a basis of 12 orbitals. 

An example of the electron density, computed for full polarization, is reported in Fig. [7]. 

The distance of the 2DEG from the surface of the heterostructure plays an important 
role and significantly affects the cell-to-cell response function: the closer it lies to the surface, 
the higher, for a given value of the bias voltage, are the potential barriers separating the 
dots in a cell, and, as a result, the steeper and more abrupt is the response function. This is 
the prevailing effect, even though it is partially compensated by the screening action of the 
surface, which reduces the cell-to-cell interaction and increases with decreasing distance. 

In Fig. |K a) we show a polarization curve which has been obtained for the same dot 
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configuration as that in Fig. || but for different values of the 2DEG depth, from 45 to 
55 nm. As expected, the polarization curve becomes smoother for increasing depth of the 
2DEG. The screening effect from the gates can be appreciated by comparing with the results 
obtained neglecting the contribution due to the images. This is seen in Fig. §(b), where we 
report the cell-to-cell response function for the previously considered gate geometry and for 
a 2DEG depth of 50 nm, computed with and without images. The image effects decrease the 
sharpness of the response function, because the dipole moment of the driver cell is screened 
by its images. 

As stated in the introduction, we have applied our Configuration-Interaction technique 
to assess the sensitivity of a two-cell structure to fabrication tolerances. First we considered 
giving one of the holes in the gate a diameter slightly different from the others: this leads to 
a variation in the potential landscape defining the cell, and, in particular, to a modification 
of the confinement energy associated with that dot. If we reduce the diameter of one of the 
holes, the confinement energy for the dot underneath will rise by a certain amount 5E. If 
5E is larger than the electrostatic splitting AEc between the two cell configurations, the 
cell will be stuck in a state with the smaller dot empty. For SE < AEc, the cell will still 
be operational, but the cell-to-cell response function will be shifted by an amount which 
depends on the ratio of 5E to AEc- The strong nonlinearity of the response function helps 
to restore the correct polarization value along a chain of cells, as long as the shift still allows 
full polarization of the driven cell for full polarization of the driver cell. 

The tolerance on the hole diameter, admissible before unrecoverable disruption of the 
operation of two coupled cells occurs, is unfortunately very small: from our calculations 
it is about one part in 10,000, for the dot sizes considered so far. Since the tolerance is 
determined by the interplay between the electrostatic splitting energy (which has an inverse 
linear dependence on the size) and the perturbation of the confinement energy (which has 
an inverse quadratic dependence on the size), it will be even tighter for smaller cells. 

In Fig. |^(a) we report the cell-to-cell response function for a gate configuration as previ- 
ously described, with the diameter of one of the dots reduced to 0.9999 times the nominal 
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value. The depth of the 2DEG is assumed to be 50 nm and the three solid curves refer to dif- 
ferent separations between cell centers. It is clear that, for a separation between the centers 
greater than 280 nm, an error of one part in 10,000 will be sufficient to unrecoverably dis- 
rupt QCA operation. For purposes of comparison, we have reported also the results (dashed 
curves) obtained neglecting the contribution from the image charges (i.e. the screening due 
to the gates and to the assumed Fermi level pinning at the semiconductor-air interface). 

We have also investigated the sensitivity to errors in the position of the gate holes. 
With a hard- wall model, such as in Ref. ||, a small shift in the position of one hole would 
not have a disrupting effect, since it would cause only a proportionately small variation of 
the electrostatic splitting and would not affect the confinement energy in any way. With a 
realistic model, as here, the situation is quite different: the confinement potential for each dot 
is determined not only by the corresponding hole, but also by the other holes belonging to the 
same cell. This means that by shifting a hole away from its nominal position, the potential 
landscape defining the dot underneath will be distorted, and the confinement energy will 
change. As a consequence, cell operation will be disrupted for unexpectedly small errors 
in hole positioning. In Fig. ^|(b) we report the cell-to-cell response function obtained for 
a 0.0275 nm shift down and to the left, of the bottom left dot of the driven cell. For the 
rest, the cell configuration is unchanged: dot diameter of 90 nm, dot-dot separation within 
a cell of 110 nm and distance between cell centers of 280 nm. Such a small displacement 
is sufficient to significantly shift the response function; a 0.05 nm shift leads already to the 
breakdown of cell operation. 

B. Multiple independent gates with adjustable voltages 

From the results shown above, it is clear that a simple hole-array implementation of 
QCA cells leads to unrealisable requirements on fabrication tolerances. This is the reason 
why we have also investigated alternative gate layouts such as that sketched in Fig. |j. If dot 
confinement is obtained via multiple independent gates, it is indeed possible to compensate 
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for geometrical tolerances by adjusting the gate voltages. 

In Fig. [XOj we show the cell-to-cell response function for a cell defined with bias voltages 
of —1.8 V applied to all gates, except for gates 2 and 6, which are fixed at -1.6 V; these bias 
values have been chosen within a reasonable voltage range so as to get four clearly confined 
dots. The separation between the centers of the driver and the driven cells is 280 nm. Also 
in this case, different depths for the 2DEG have been considered, from 40 to 55 nm. 

A cross-section of the confinement potential, cut across the two upper dots, is shown in 



Fig. [TTjfor a 2DEG at 35 nm (solid line), 45 nm (dotted line) and 55 nm (dashed line): the 
height of the barriers is relatively low, but their width guarantees a low enough transparency 
for strong localization of the electrons and thereby correct operation of the cell. 

We have also investigated the dependence of the cell-to-cell response function on electrical 
asymmetries, for the case of a 2DEG 50 nm depth and a distance of 280 nm between 
the centers of the driver and the driven cell. For this purpose, the voltage applied to 
gate 3 has been made slightly more negative by an amount SV. The results for 5V = 
—0.05, —0.1, —0.2 mV are shown in Fig. |T2"|(a). The most visible effect is a shift in the cell- 
to-cell response function, which is somewhat proportional to the variation in the applied 
voltage. This does not disrupt the operation of a QCA chain, as long as full polarization of 
one cell can produce full polarization of the neighboring cell. Therefore, for this particular 
cell, we expect a maximum tolerance on the gate voltages of about 0.4 mV. This may seem 
difficult to achieve at first sight, but it is important to keep in mind that it is a shift between 
gate biases: larger variations in the overall average value of the gate voltages are allowed, 
as long as they do not alter cell occupancy. 

As previously mentioned, the screening effect due to the gates and to the charge at 
the semiconductor-air interface decreases the strength of the electrostatic interaction and, 
therefore, the energy splitting between the two possible cell polarizations. Hence, the effects 
of asymmetry decrease with increasing depth of the 2DEG, as the image effects are reduced. 
This phenomenon is clearly visible in Fig. [T^(b), where the cell-to-cell response function is 
plotted for a voltage shift of —0.1 mV on gate 3 and various values of the 2DEG depth 
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ranging from 35 nm to 50 nm. The shift in the response function decreases as the depth of 
the 2DEG increases, and the effect of the image term becomes less important. 

Fabrication tolerances would also disrupt the operation of this type of cell, but they can 
be compensated for by fine adjustments to the bias voltages applied at the gates. As an 
example, we consider a cell with gate 3 shifted by 5 nm to the right. An iterative procedure 
was developed for computing the new bias voltages that will restore the symmetry of the 
structure. For a perfectly symmetric structure there is a four-fold quasi-degeneracy of the 
one-electron states, corresponding to the four-fold symmetry of the cell. When symmetry is 
disrupted, this quasi-degeneracy is lifted, and the first four eigenvalues differ from each other 
by an amount which is no longer negligible. Our strategy is to evaluate the difference between 
the first and the fourth eigenvalue, and then adjust each gate voltage in such a way as to 
minimize this difference. We vary one gate voltage at a time, from gate 1 to gate 7 (with the 
exception of gate 4, which does not affect cell symmetry), and then repeat the cycle, starting 
again from gate 1, until the splitting between the first and the fourth eigenvalue is smaller 
than an assigned threshold. In order to avoid getting stuck in local minima, it is convenient 
to perform the whole minimization procedure several times, for values of the displacement 
of gate 3 increasing with a geometric progression from 0.01 nm to 5 nm. Finally, small 
adjustments are made manually, until a symmetric cell-to-cell response function is obtained. 
The gate voltages needed to symmetrize the cell are listed in Table |, while the cell response 
function for the symmetrized cell is shown in Fig. [H| This demonstrates that a 10% error 
in the position of one of the gates can be fully compensated. State of the art fabrication 
techniques allow geometrical tolerances of this order of magnitude or smaller, thus such a cell 
is actually manuf act ur able, although not useful for large scale applications where it would be 
impossible to tune each cell separately. By acting on the gate voltages it is also possible to 
compensate for the presence of randomly distributed stray charges, which would also disrupt 



QCA operation. We notice that the cell-to-cell response function of Fig. [13] is steeper than 
that for a geometrically symmetric cell with analogous parameters (see Fig. |Tt]). This is 
a consequence of a slight variation in the barrier heights and widths (due to the different 
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applied voltages) and of the exponential dependence of the tunneling coefficients between 
the dots on such barrier characteristics. A chain of such cells can be fabricated by repeating 
this same gate layout in the horizontal direction. However, lateral branching, needed for 
the implementation of logic gates U, is not allowed, due to the lateral extension of the 
leads required for feeding the bias voltages. From this point of view, a more promising 



implementation would be that suggested by Chen and Porod |L2| , with central enhancement 
gates in each dot: by adjusting the voltage of such a central gate it would be possible to 
correct asymmetries, while keeping the possibility of lateral branching. This implementation, 
however, poses serious fabrication problems, because of the difficulties involved in separately 
contacting all of the central gates. 

We emphasize that the compensation procedure described above is not proposed as a 
practical method, rather our aim has been to demonstrate that compensation is possible in 
principle. 



V. CONCLUSIONS 

We have shown that the Configuration Interaction method allows us to solve the many- 
electron Schrodinger equation for a QCA cell made up of four coupled dots. This represents 
a definite advance over the state of the art in the simulation of multiple quantum dot 
systems: it is possible to include realistic confinement potentials without resorting to some 
type of mean-field approximation (local density approximation, Hartree-Fock, Hartree-Fock- 
Roothaan), which fail to converge when single-electron states are strongly degenerate and the 
electrostatic interaction is comparable to the confinement energy. A Hubbard-like approach 
to QCA cells is also feasible and has been shown to provide a qualitative understanding of 
the underlying physics ||, but requires a set of phenomenological parameters, such as the 
on-site electrostatic interaction, the dot confinement energy and the tunneling energy, which 
cannot be easily obtained from the geometrical structure and from experiments. With the 
CI method it is sufficient to determine the confinement potential from the layer structure 
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and the gate layout. This is an important advantage, which makes simulations based on the 
CI method a reliable and effective design tool. 

A disadvantage of the CI approach lies in the large computational resources which are 
required, since the number of Slater determinants to be considered exhibits a combinatorial 
increase with the number of electrons in the system. In this paper, we have presented results 
for two electrons, and a QCA cell with up to six electrons is currently being investigated. 
Calculations for a larger number of electrons would require prohibitive memory sizes (well 
above 1 Gbyte) or extremely long computation times. 

We have focused our investigation on the sensitivity to fabrication tolerances for two 
coupled QCA cells. Our results demonstrate that the implementation of a simple "hole- 
array" approach is not feasible, because it would require a precision in the diameter of each 
hole that is well beyond the present state of the art in electron-beam lithography. 

We have proposed an alternative cell layout, based on seven gates, whose bias voltages 
can be independently adjusted: this approach is within the capabilities of current fabrication 
technologies, since geometrical errors in the gate positions can be corrected by means of ap- 
propriate voltage variations. Admissible voltage tolerances are rather small, but achievable 
using resistive voltage dividers cooled down together with the sample. The approach we are 
proposing should allow an experimental demonstration of the QCA principle, from a single 
cell up to a chain of cells, but it is not suitable for the realization of logic gates, due to 
the impossibility of lateral branching, which is prevented by the leads reaching each gate. 
It is also impractical for large-scale integration, due to the need for individual adjustment 
of each single cell. Viable logic circuits will require drastically different solutions and new 
architectural concepts. 
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APPENDIX A: 



The problem of computing the matrix elements of the Hamiltonian in Eq. (H) between 



two Slater determinants is well known [19]. For the diagonal elements one finds 



i 

+ \ Y,((PnikVn jk \gWn ik tPn jk ) ~ {<Pn ik <Pn jk \g\<Pn jk <fn ik )) , (Al) 

y 

where in general 

(<Pi<Pj\g\<Pi<p m ) = J dqi dq 2 wl{qi)tf{q2)g{rur 2 )(pi{qi)ip m {(l2) • (A2) 

As far as the computation of the off-diagonal matrix elements of the Hamiltonian of 
Eq. (^|) between two different Slater determinants $fc , is concerned, there are some 
"selection rules" (Slater's rules 0]) which state that there are only two possible cases in 
which (<&k\H\<&w) is not vanishing, i.e. when , either differ by one single spin-orbital 
or by two: 

1. one spin-orbital difference (<f nik 7^ fn-u) 

{$ k \H\$ v ) = ((p nik \h\cp n . k> ) + 

^2{( ( Pr H k ( Pn j j,\g\(p nik ,(P njk ) ~ (Vn ik <Pn jk \g\<Pn jk <Pn ik ,}) , (A3) 

2. two spin-orbital difference (^„ ife ^ ^> n%k , and (p njk ^ ip n . y ) 

(Q k \H\$ w ) = (cp n . k ip n . k \g\ip n . kl {p n ) - ((p nik ip n . k \g\(p n cp n . k ,} . (A4) 



The expressions in Eqs. ( |A3| , A4|) refer to the case in which the spin-orbitals that are 



common to both Slater determinants occur in the same columns. If this is not the case, it 
is possible to perform a permutation of the columns of one determinant, so that the above 
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condition is satisfied; the permutation has the effect of changing the sign of the matrix 
element if it is an odd order permutation. 

Finally, it is worth mentioning that Eqs. (|A1| , |A3| , |A4j ) are valid only if the orthonormality 
condition on the spin-orbitals is satisfied. 
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FIGURES 

FIG. 1. Schematic representation of two coupled QCA cells: tunneling of electrons is possible 
along the dashed lines. 

FIG. 2. Subdivision of a cell into four quadrants. 

FIG. 3. QCA cell obtained by depleting a two-dimensional electron gas by means of a metallic 
gate with four holes, which define the four quantum dots. 

FIG. 4. Gate layout for the definition of four independently adjustable quantum dots; all 
distances are in nanometers. 

FIG. 5. Confinement potential at a depth of 50 nm from the surface of the heterostructure, 
produced by a gate with four holes of diameter 90 nm and placed at the corners of a 110 nm square. 

FIG. 6. Cell-to-cell response function obtained for a cell defined by a gate, kept at —0.5 V, 
with four 90 nm holes placed at the corners of a 110 nm square. The depth of the 2DEG is 50 nm 
and the separation between the centers of the driven and the driver cells is 400 nm. The different 
symbols indicate different numbers of basis functions. 

FIG. 7. Electron density, in arbitrary units, for a completely polarized cell defined by a gate, 
kept at —0.5 V, with four 90 nm holes placed at the corners of a 110 nm square. The depth of the 
2DEG is 50 nm. 

FIG. 8. Cell-to-cell response functions obtained for a cell defined by a gate, kept at —0.5 V, 
with four 90 nm holes placed at the corners of a 110 nm square. The separation between the centers 
of the two cells is 400 nm. (a) Comparison between the response functions for different values of 
the 2DEG depth, (b) Comparison between the response functions obtained with and without the 
contribution of the image charges. 
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FIG. 9. Cell-to-cell response functions obtained for a cell denned by a gate, kept at —0.5 V, 
with four holes placed at the corners of a 110 nm square, (a) The bottom left hole has a diameter 
of 89.991 nm and the other 3 holes of 90 nm; the depth of the 2DEG is 50 nm. Results for three 
different values of the separation between the cell centers are reported, including the effect of image 
charges (solid lines) or neglecting it (dashed lines), (b) The four holes have a diameter of 90 nm, 
but the bottom left hole has been moved down and to the left by 0.0275 nm. 

FIG. 10. Cell-to-cell response functions obtained for a cell defined with the seven-gate layout 
and a distance between cell centers of 280 nm. The gates are kept at —1.8 V, except for gates 2 
and 6, which are kept at —1.6 V. Results for a depth of the 2DEG varying between 40 and 55 nm 
have been reported. 

FIG. 11. Cross-section of the confining potential obtained with the seven-gate layout, cut across 
the two upper dots, at a depth from the surface of 35 nm (solid line), 45 nm (dotted line), and 
55 nm (dashed line). 

FIG. 12. Cell-to-cell response functions obtained for a cell defined with the seven-gate layout 
and for a distance between cell centers of 280 nm. The gates are kept at —1.8 V, except for gates 
2 and 6, which are kept at —1.6 V, and gate 3, which is kept at —1.8 V+5V. (a) The depth of the 
2DEG is 50 nm and the results for three different values of SV have been reported, (b) The value 
of SV is —0.1 mV and results for different depths of the 2DEG are shown. 

FIG. 13. Cell-to-cell response functions obtained for a cell defined with the seven-gate layout, 
a distance between cell centers of 280 nm, and a 2DEG depth of 50 nm. Gate 3 has been shifted to 
the right by 5 nm, and the bias voltages are those listed in Table 1, chosen to restore cell symmetry. 
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TABLES 

TABLE I. Values of the bias voltages to be applied to the gates defining the geometrically 
asymmetric cell shown in Fig. 4, in order to symmetrize it. 

Gate Bias voltage 

Vi -1.665577 V 

V 2 -1.728683 V 

V 3 -1.845800 V 

Vi -1.800000 V 

V 5 -1.807502 V 

V 6 -1.592665 V 

V 7 -1.793715 V 
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